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ORIGINAL ARTICLE 

Follow-up of loci from the International Genomics of 
Alzheimer's Disease Project identifies TRIP 4 as a novel 
susceptibility gene 

A Ruiz^'^^ S Heilmann^'^'^^ T Becker'^'^'^^ I Hernandez\ H Wagner^ M Thelen^ A IVlaule6n\ M Rosende-Roca\ C Bellenguez^'^'^ 
JC Bis^° D Harold", A Gerrish", R Sims", O Sotolongo-Grau\ A Espinosa\ M Alegret\ JL Arrieta^^ A Lacour^ M Leber^ J Becker^ 
A Lafuente\ S Ruiz\ L Vargas\ O Rodriguez\ G Ortega\ M-A Dominguez\ IGAP^^ R Mayeux^^'^^ JL Haines^^'^^ MA Pericak-Vance^^'^^ 
LA Farrer^^'2°'2^'22'2^ GD Schellenberg^^ V Chouraki^^ LJ Launer^^ C van Duijn^^'^^'^^ S Seshadri^^ C Antunez^^ MIVl Breteler^ 
M Serrano-Rios^°, F Jessen'^'^ L Tarraga\ MIVl Nothen^'^ W Maier'^'^ M Boada^'^^ and A Ramirez^'^ 

To follow-up loci discovered by the International Genomics of Alzheimer's Disease Project, we attempted independent replication 
of 19 single nucleotide polymorphisms (SNPs) in a large Spanish sample (Fundacio ACE data set; 1808 patients and 2564 controls). 
Our results corroborate association with four SNPs located in the genes INPP5D, MEF2C, ZCWPWl and FERMT2, respectively. Of 
these, ZCWPWl was the only SNP to withstand correction for multiple testing (P = 0.000655). Furthermore, we identify TRIP4 
(rs74615166) as a novel genome-wide significant locus for Alzheimer's disease risk (odds ratio = 1.31; confidence interval 95% 
(1.19-1.44); P = 9.74x10"^). 
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INTRODUCTION 

Alzheimer's disease (AD) is a complex multifactorial neuropsychia- 
tric disorder whose etiology involves both environmental and 
genetic factors. The major genetic risk factor for AD is the 
apolipoprotein E4 {AP0E4) allele.^ For 17 years after its discovery, 
this remained the only confirmed genetic risk factor for the 
disorder. However, subsequent meta-analyses of genome-wide 
association studies identified further genetic risk loci. These 
include signals close to, or within, candidate genes such as 
CLU,^'^ PICALM,^ CRh^ BlNlf^ ABCA7,^ EPHAh^'^ CD33,^'^ CD2AP^'^ 
and ATP5H/KCTD2/ as well as the MS4A gene cluster.^'^'^ 

The International Genomics of Alzheimer's Disease Project 
(IGAP) is the largest genetic epidemiology investigation of AD risk 



to date. In 2013, the IGAP reported a mega meta-analysis, which 
was divided into a discovery step (stage 1) and a replication step 
(stage 2). This mega meta-analysis comprised 74 046 samples, 
including those of 25 580 AD cases, and identified 1 1 new loci, 
thus doubling the number of genome-wide significant loci 
reported for AD.^ The analysis also identified 13 suggestive loci. 
These findings may serve as the starting point for novel 
discoveries in future AD genomics studies. 

Four of the 1 1 genome-wide significant loci in the IGAP analyses 
reached significance in stage 1 (rs8093731 DSG2; rs28834970 
PTK2B; rsl 1218343 SORLl; rs10498633 SLC24A4). The remaining 
seven only reached genome-wide significance in stage 2 — that is, 
after the inclusion of the replication sample (rs35349669 INPP5D; 
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Table 1. Results for the 19 IGAP SNPs in the Fundacio ACE data set 
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rs190982 MEF2C, rs2718058 A/MES; rs1476679 ZCWPWl; MATERIALS AND METHODS 

rsl 0838725 CELFV, rsl 71 25944 FERMT2, and rs7274581 CASS4). Patients and controls 

Replication of the IGAP findings in an independent series is The present study involved 4372 individuals. These included 1 808 patients 

therefore warranted. The present report describes the follow-up of ^ith a possible or probable diagnosis of AD, as assigned by a 

10 novel and nine suggestive IGAP loci using subjects drawn from neurologist,^'^ ° and 2564 unrelated healthy controls from the Spanish 

the Spanish Fundacio ACE cohort.^ general population who were selected from the Neocodex bio-bank." The 
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AD cases were recruited consecutively from three centers: Barcelona 
(n = 1 627); Madrid (n = 1 61 ); and Murcia [n = 20). None of these AD patients 
had been included in the IGAP replication analyses. To avoid the issue of 
population stratification, all cases and controls were of Spanish ancestry, 
which was defined as a history of two generations of registered Spanish 
ancestors. The demographic characteristics of the Fundacio ACE partici- 
pants are described elsewhere.^ Written informed consent was obtained 
from all participants, or from their legal representatives when necessary. 
The study was approved by the respective ethics committees, and was 
performed in accordance with the Code of Ethical Principles for Medical 
Research Involving Human Subjects of the World Medical Association. 

DNA isolation and genotyping 

DNA was extracted using 'Magnapure' technology (Roche Diagnostics, 
Mannheim, Germany). Twenty-five single nucleotide length polymorph- 
isms (SNPs) with an AD-association risk of below P< 10"^ in the IGAP 
consortium study were selected for replication. The primer molecules for 
the multiplex reaction were designed using the Assay Design Suite tool 
(www.mysequenom.com, Sequenom, San Diego, CA, USA). Assay designs 
were successful for 21 of the 25 selected variants. Four SNPs (rs72807343 
(SQSTM1); rs9271192 (HLA-DRB5/HLA-DRB1); rs2337406 (IGH@); and 
chrl 7:61,538,148 (ACE)) were rejected during this phase due to technical 
problems. Primer sequences and assay conditions for the genotyped SNPs 
are available upon request. 

Quality control 

A total of 1808 AD patients and 2564 controls were genotyped for 21 SNPs 
using Sequenom's Mass Array System (Sequenom) and iPlex Gold reagents 
in accordance with the manufacturer's instructions. Only SNPs with a call 
rate of ^95% and a Hardy-Weinberg equilibrium P-value of >0.01 in the 
whole data set were included in the subsequent analyses (Supplementary 
Table 1 ). All SNP major and minor alleles and allelic frequencies obtained in 
the Fundacio ACE data set were fully consistent with those reported by the 
IGAP consortium (Table 1). The overall conversion rate was 96.7%. The 
SNPs rsl 0751 667 {AP2A2) and rsl 0838725 {CELFl) failed quality control and 
were excluded from the statistical analyses. The 19 successfully genotyped 
SNPs and their status in the IGAP analyses (that is, genome-wide 
association studies significant or suggestive) are specified in Table 1. 

Statistical analysis 

Genetic association analyses and calculation of allelic frequencies and 
Hardy-Weinberg equilibrium were conducted using the online tool at the 
TUM Helmholtz Center (Munich, Germany; http://ihg.gsf.de/cgi-bin/hw/ 
hwa1.pl). Briefly, tests adapted from those of Sasieni^^ were used. Age- and 
sex-adjusted binary logistic regression analyses were performed using SPSS 
15.0 software (SPSS, Chicago, IL, USA). In addition, Mantel-Haenzsel- 
stratified analyses were conducted according to gender and the presence 
or absence of the AP0E-e4 allele (Supplementary Tables 2 and 3, 
respectively). Breslow-day tests were conducted to measure the signifi- 
cance of SHP xAPOE and SNP x gender interactions. All SNP calculations 
were double-checked using PLINK or INTERSNP software.^ ^'^"^ Meta-analyses 
were conducted using the PLINK software. All results were doubled- 
checked using Ken Rothman's Episheet spreadsheet and PLINK (for details 
see http://pngu.mgh.harvard.edu/~purcell/plink/; http://krothman.hostbyet2. 
com/episheet.xis). Power calculations were performed using Episheet. 



RESULTS AND DISCUSSION 

In the present replication effort, a nominally significant signal 
(P<0.05) was detected for six of the 19 investigated SNPs: 
rs35349669 at INPP5D (P = 0.023); rsl 90982 at MEF2C (P= 0.0062); 
rsl 476679 at ZCWPWl (P = 0.00065); rsl 71 25944 at FERMT2 
(P = 0.018); rs8035452 at SPPL2A (P- 0.031); and rs74615166 at 
TRIP4 (P = 0.0032) (Table 1). Of these, rs1 476679 at the ZCWPWl 
locus, which had shown genome-wide significance in stage 2 of 
the IGAP analyses, was the only SNP to withstand correction for 
multiple testing (P = 0.000655). However, the observed inflation 
factor for was X = 4.7. As in genome-wide analyses, the X-value 
was computed as the median of the x^-test statistics obtained for 
the 19 investigated SNPs. Under the null hypothesis of no 
association, the expected A value is 1 in the absence of true 



association. In our study, however, the A value is 4.7, which is a 
strong sign of an overall increased degree of significant 
associations for the SNPs investigated in the Spanish cohort. In 
the present context, the A value is considered to indicate the 
overall degree of positive associations for a small set of SNPs, in 
contrast to the genome-wide setting, where it is used as an 
indicator of residual inflation caused by spurious association. Thus, 
the present observations in our study are unlikely to represent 
chance findings. Furthermore, four of these five nominally 
significant association signals displayed the same effect direction 
as that reported by the IGAP (Table 1). The exception was the 
marker rs8035452. The IGAP reported this as a suggestive signal. 
However, an effect in the opposite direction was found in the 
Fundacio ACE data set. This observation might reflect a lack of 
power in our data set to detect this signal. Alternatively, the 
original finding may represent a false-positive. 

No significant association was found for three of the four 
genome-wide significant loci detected during IGAP Stage 1 
(rsl 1218343 at SORLl; rs10498633 at 5LC24A4; rs8093731 at D5G2). 
However, since the effect sizes and directions of these three loci 
were fully consistent with those reported by the IGAP, our failure to 
replicate them may have been attributable to a lack of power. 

In total, seven of the 10 investigated genome-wide significant 
loci from the IGAP displayed a consistent effect in the present data 
set. The non-consistent effects observed for the NME8, PTKB2, and 
CASS4 signals may have been attributable to a lack of power. It 
should be noted that, although our series may appear under- 
powered compared with the IGAP data set, the present study had 
on average a power of 45% to detect each of the genome-wide 
association studies significant signals reported by the IGAP. 

The results of the APOE and gender-adjusted stratified analyses 
suggested that for most of the 19 investigated SNPs, APOE status 
and gender had little impact on effect size or the association 
results (Supplementary Table 3). Interestingly, nominal P-values for 
APOE interaction were obtained for two nonsignificant SNPs in our 
series (rs7295246 ADAMST20; and rs7225151 SCIMP). Both loci 
were reported as being suggestive by the IGAP. The results of the 
present stratification analyses support the hypothesis that these 
two loci represent susceptibility factors in only a fraction of AD 
patients, and that their effects are dependent upon the AP0E-E4 
genotype. This observation may facilitate determination of their 
role in AD development in future studies (Supplementary Table 3). 

Of the nine suggestive loci proposed by the IGAP, only one SNP 
was significant in the present analyses (Table 1). The statistically 
significant signal was obtained for rs74615166 at the thyroid 
receptor interacting protein gene 4, TRIP4, locus (odd ratio = 1.519 
(1.148-2.012), P = 0.0032). This variant had a minor allele 
frequency of 0.02 in both the IGAP and the Fundacio ACE 
(TalDle 1). Interestingly, a larger effect size was observed in the 
Fundacio ACE data set than in the IGAP. However, an advantage of 
the present analyses was that this SNP was genotyped directly, 
whereas the IGAP had to rely in part on imputed genotypes. Since 
imputation for rarer variants is more difficult, this might explain 
why a stronger effect was observed in the present cohort. The 
present findings for the nine suggestive IGAP signals may indicate 
that these loci have a weaker effect on AD risk than the genome- 
wide significant SNPs. As a direct consequence, the power to 
detect them using our data set is relatively low (33% on average 
for suggestive signals). However, our data set had a >99.9% power 
to detect at least one suggestive locus (0.33),^ thus explaining the 
results for the TRIP4 locus. 

Meta-analysis of the present results with the IGAP meta-analysis 
data identified TRIP4 as a novel genome-wide significant locus. 
The new susceptibility AD SNP is located within the eleventh 
intron of TRIP4 (15q22.31; rs74615166; OR =1.31 (1.17-1.42), 
P = 9.74x 10~^; Table 1). According to publicly available databases 
(genome.cse.ucsc.edu), TRIP4 is highly expressed in the immune 
system and has been detected in various tissues, including the 
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brain. Research suggests that TRIP4 is a component of the nuclear 
receptor-coupled co-activation machinery that enables or disables 
DNA transcription.^^ A homolog of the TRIP4 gene in Caenorhab- 
ditis elegans showed elevated transcript levels in aged or starved 
adults, which suggests that TRIP4 has a role in cellular 
maintenance or survival. The TRIP class of proteins show thyroid 
hormone-dependent interaction with their receptors, and the 
association between TRIP4 and AD risk may partly explain previous 
findings of an association between low thyroid-stimulating 
hormone levels in clinically euthyroid subjects and increased AD 
risk.^^ TRIP proteins show a similar ligand-dependent interaction 
with the retinoid X receptor. This is of interest, since a recent AD 
mouse model study reported that administration of the retinoid X 
receptor agonist bexarotene resulted, within hours, in enhanced 
clearance of soluble A|3.^^ Besides TRIP4, the linkage disequili- 
brium block that contains rs74615166 includes several other 
candidate genes, such as CSNKIGI. CSNKIGI is a member of the 
CK-1 family, and its gene product has been implicated in the 
amyloid cascade.^^ 

A major limitation of the present study was the lack of power to 
confirm all true associations. Therefore, our negative results 
cannot be interpreted as confirmation of a lack of association 
for the respective SNPs, which remain putative susceptibility loci 
for AD. A fraction of the suggestive SNPs reported by the IGAP 
may be genuine, and these SNPs warrant further investigation. 

The results obtained for TRIP4 underscore the importance of 
follow-up and comprehensive replication of consortia results. 
Further genotyping and re-sequencing efforts to investigate TRIP4 
and the other IGAP loci are underway in order to elucidate the role 
of TRIP4 in AD risk and corroborate further genuine signals. 
Further studies are now warranted to identify the functional 
mechanism underlying the association between TRIP4 and AD. 
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